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Abstract. In this paper we give an algorithmic characterization of rank two 
locally nilpotent derivations in dimension three. Together with an algorithm 
for computing the plinth ideal, this gives a method for computing the rank of 
a locally nilpotent derivation in dimension three. 



1. Introduction 

Let JC be a commutative field of characteristic zero, JC>- n ' be the ring of poly- 
nomials in n variables with coefficients in /C and AutjcQO- n ') be the group of K- 
automorphisms of /C'™'. The structure of Aut^KP^) is well understood [21]. How- 
ever, remains a big mystery for n > 3. 

In order to understand the nature of Autjc(IC^) it is natural to investigate 
algebraic group actions on the affinc n-space over K. Actions of the algebraic 
group (JC,+) on affine spaces are commonly called algebraic G a -actions, and they 
are all of the form exp(tX) t £ic where A" is a locally nilpotent ^-derivation of the 
polynomial ring £M . 

Locally nilpotent derivations of /C' 2 ' are completely classified, and this classifi- 
cation is algorithmic [25]. In dimension three, D. Daigle, G. Freudenburg and S. 
Kaliman obtained several deep results which constitute a big step towards a classi- 
fication, see [17] and the references therein. However, some of these results which 
are obtained by using topological methods are not of algorithmic nature. It would 
of course be very nice to obtain an algorithmic classification of locally nilpotent 
derivations in dimension three, but this seems to be a difficult problem. This paper 
addresses the less ambitious problem of computing some invariants, namely the 
plinth ideal and the rank, of a locally nilpotent derivation in dimension three. 

The paper is structured as follows. In section 2 we set up notation and give 
the main results to be used. Section 3 concerns minimal local slices and how to 
compute them in dimension 3. This gives an algorithm for computing a generator 
of the plinth ideal of a locally nilpotent derivation in dimension three. In section 
4 we give an algorithm to compute the rank of a locally nilpotent /C-derivation in 
dimension 3. 
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2. Notation and basic facts 

Throughout this paper /C is a commutative field of characteristic zero, all the 
considered rings are commutative of characteristic zero with unit and all the consid- 
ered derivations are nonzero. A derivation of a /C-algebra A is called a ^-derivation 
if it satisfies X(a) = for any a G /C. 

2.1. Classical results on locally nilpotent derivations. A derivation X of a 
ring A is called locally nilpotent if for any a G A there exists a positive integer n 
such that X n (a) = 0. An element s of A satisfying X(s) ^ and X 2 (s) = is 
called a local slice of X. If moreover X(s) = 1 then s is called a slice of ,Y. A 
nonzero locally nilpotent derivation needs not to have a slice but always has a local 
slice. 

The following result, which dates back at least to [27], concerns locally nilpotent 
derivations having a slice. 

Lemma 2.1. Let A be a ring containing Q and X be a locally nilpotent derivation 
of A having a slice s. Then A = A x [s] and X = d s . 

Let A be a ring and X be a derivation of A. The subset {a G A; X(a) = 0} of A 
is in fact a subring called the ring of constants of A" and is denoted by A x . When 
A is a domain and X is locally nilpotent, the ring of constants A x is factorially 
closed in A, i.e., if a G .A* and a = 6c then b, c G A 1 *. Also, the fact that A is of 
characteristic zero implies that A* C A x '. 

Locally nilpotent derivations in two variables over fields of characteristic are 
well understood. We have in particular the following version of Rentschlcr's theorem 
[25]. 

Theorem 2.2. Let X be a locally nilpotent IC-derivation oflC[x, y\. Then there exist 
two polynomials f,g of JC[x,y] and a univariate polynomial h such that JC[f,g] = 
K[x,y], )C[x,y] x = £[/] and X = h(f)d g . 

As a consequence of theorem 2.2, if A is a UFD containing Q and X is a locally 
nilpotent A-dcrivation of „4[x,?/] then there exists / G A[a;,y] and a univariate 
polynomial h such that A[x, y] x = A[f] and X = h(f)(d y fd x — d x fd y ), see [9, 12]. 

In the case of /C^ we have the following result proved by Miyanishi [24] for the 
case K, = C and may be extended to the general case in a straightforward way by 
using Kambayashi's transfer principle [20], see also [6] for an algebraic proof. 

Theorem 2.3. Let X be a locally nilpotent IC-derivation of K[x,y,z\. Then there 
exist polynomials f,g such that JC[x , y , z] x = JC[f,g]. 

Contrary to lemma 2.1 and theorem 2.2 which are of algorithmic nature, it is not 
clear from the existing proofs of theorem 2.3 how to compute, for a given locally 
nilpotent derivation X of )C{x,y,z], two polynomials /, g such that K.[x, y, z] x = 
K\f,g\. 

2.2. Coordinates. A polynomial / G K\xi, . . . ,x n ] is called a coordinate if there 
exists a list of polynomials f±,..., f n -i such that K\x\, . . . , x n ] = IC[f, fi, . . . , f n -i]- 
In the same way a list /i, . . . , f r of polynomials, with r < n, is called a system of co- 
ordinates if there exists a list f r +i, ■ ■ ■ , fn of polynomials such that K.[xi, . . . , x n ] = 
/C[/i, . . . , f n ]- A system of coordinates of length n is called a coordinate system of 
K[xi, ...,x n ]. 
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The famous Abhyankar-Moh theorem [1] states that a polynomial / in K\x, y) is 
a coordinate if and only if JC[x,y]/f is /C-isomorphic to KP^. In the case of three 
variables we have the following result proved by Kaliman in [19] for the case K. = C 
and extended to the general case in [10] by using Kambayashi's transfer principle 
[20]. 

Theorem 2.4. Let f be a polynomial in K.[x,y,z] and assume that for all but 
finitely many a £ K. the K-algebra K[x, y, z]/(f — a) is IC-isomorphic to JO 2 ' . Then 
f is a coordinate of JC[x, y, z] . 

As a consequence of theorem 2.4, if a polynomial / is such that K.(f)[x,y, z] 
is /C(/)-isomorphic to )C(f)^ then / is a coordinate [23, 13]. This is the version 
of Kaliman's theorem we will use in this paper. However, it is not clear how to 
compute polynomials g, h such that /, g, h is a coordinate system of /C[x, y, z] since 
the original proof given in [19] is of topological nature. 

2.3. Rank of a derivation. Let X be a /C-derivation of K\x\ = K\x\, ■ ■ • , x n ]. As 
defined in [16] the co-rank of X, denoted by corank(X) , is the unique nonnegativc 
integer r such that /C[x]' t ' contains a system of coordinates of length r and no system 
of coordinates of length greater than r. The rank of X, denoted by rank(X), is 
defined by rank(X) = n — corank(X). Intuitively, the rank of X is the minimal 
number of partial derivatives needed to express X. The only one derivation of rank 
is the zero derivation. Any /C-derivation of rank 1 is of the form p(fi, . . . , f n -i)df n , 
where f\ , . . . , f n is a coordinate system. Such a derivation is locally nilpotent if 
and only if p does not depend on /„. 

Let X be a nonzero locally nilpotent /C-derivation of K[x\, . . . ,x n ], write X = 
aid Xi and set c = gcd(ai, . . . , a n ). We say that X is irreducible if c is a constant 
of K,* . It is well known that X(c) — and X = cy where y is an irreducible 
locally nilpotent /C-derivation. Moreover, this decomposition is unique up to a 
unit, i.e., if X = Ci^i, where is irreducible, then there exists an nonzero constant 
/.t G K* such that C\ = /ic and y = /iJV Given any irreducible locally nilpotent 
/C-derivation of K\x\, . . . , x n ] and any c such that X(c) = 0, the derivations X and 
cX have the same rank. Thus, for rank computation we may reduce, without loss 
of generality, to irreducible derivations. 

2.4. The plinth ideal of a locally nilpotent derivation. Let A be a ring, X 
be a locally nilpotent derivation of A and let 

S x := {X(a) | X 2 (a) = 0}. 

Clearly S x is an ideal of A x , called the plinth ideal of X. In case A = K.[x,y, z] 
we have the following result from [10] which is a consequence of faithful flatness of 
JC[x,y, z] as K.[x,y, z] x -module. 

Theorem 2.5. Let X be a locally nilpotent IC-derivation of IC[x,y, z]. Then the 
ideal S x is principal. 

As we will see in the sequel, the ideal S x contains a crucial information for 
computing the rank of a locally nilpotent derivation in dimension 3. Before going 
on to the details on how to exploit this information for our purpose we first focus 
on how to compute a generator of this ideal. 



4 



M. A. BARKATOU, H. EL HOUARI, AND M. EL KAHOUI 



3. Minimal local slices 

In this section we give an algorithm to compute a generator of the ideal S x in 
case X is a locally nilpotent /C-derivation of JC[x,y, z\. Besides theorem 2.5, our 
algorithm strongly depends on the fact that /C[a;, y, z] x is generated by 2 polynomi- 
als. Since we do not have at disposal an algorithmic version of Miyanishi theorem 
we assume a generating system of K[x, y, z] to be available. 

Definition 3.1. Let A be a domain and X be a locally nilpotent derivation of A. 
A local slice s of X is called minimal if for any local slice v such that X(v) \ X{s) 
we have X(v) = /iX(s), where /i is a unit of A. 

Lemma 3.2. Let A be a domain and X be a locally nilpotent derivation of A. Let 
s be a local slice of X , p be a factor of X(s) = c and write c — pc^. Then there 
exists si G A such that X(si) = C\ if and only if the ideal pA contains an element 
of the form s + a where a G A . 

Proof. "=>") Assume that there exists a local slice si of X such that X(s\) = c\. 
Then X{ps\ — s) = and so ps\ — s = a where a is a constant of X. This proves 
that pA contains s + a. 

) Assume now that the ideal pA contains an element of the form s + a, 
where a is a constant of X, and write s + a = ps\. Then X(s) = pX(si) and so 
X(s 1 )^c 1 . □ 

Proposition 3.3. Let A be a UFD, X be a locally nilpotent derivation of A and s 
be a local slice of X . Then the following hold: 

i) there exists a minimal local slice sq of X such that X(so) \ X(s), 

ii) in case S x is a principal ideal, it is generated by X(s) for any minimal local 
slice s of X . 

Proof, i) Let ,s be a local slice of X and write X(s) = fip" 11 ■ ■ -p™ r , where /i is 
a unit and the p^s arc primes, and set m = J^- m^. We will prove the result by 
induction on m. 

For m = we have X(s) = fi, and so /^ _1 s is a slice of X. This shows that s is 
a minimal local slice of X. Let us now assume the result to hold for m — 1 and let 
s be a local slice of X, with X(s) = /ip™ 1 • • -p™ T and ^ rrii = m. Then we have 
one of the following cases. 

- For any i = 1, . . . ,r the ideal piA does not contain any element of the form 
s + a with X(a) = 0. In this case s is a minimal local slice of X by lemma 3.2. 

- There exists i such that PiA contains an element of the form s + a, with 
X(a) = 0. Without loss of generality we may assume that i = 1. If we write 
s + a = p\S\ then X(si) = p™ ll ~ p™ 2 ■ ■ ■p™ r , and by using induction hypothesis we 
get a minimal local slice sq of X such that X(sq) \ X(si). Since X(s{) \ X(s) we 
get the result in this case. 

ii) Assume now that S x is principal and let c be a generator of this ideal, with 
c = X(so). Let s be a minimal local slice of X. Since X(s) G S x we may write 
X(s) = fiX(so). The fact that s is minimal implies that [i is a unit of A x , and so 
X(s) generates S x . □ 

The main question to be addressed, if we want to have an algorithmic version of 
proposition 3.3, is to check, for a given prime p of A, whether pA Pi A x [s] contains 
a monic polynomial of degree 1 with respect to s. In case A is an affine ring over 
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a computable field K, this problem may be solved by using Grobner bases theory, 
see e.g., [2, 3, 5]. We only deal here with the case where A and A x arc polynomial 
rings over a field since this fits our need. 

Proposition 3.4. LetX be an ideal o//C[xi, . . . ,x n ] = lC[x\ and h = hi, . . .,ht be a 
list of algebraically independent polynomials of K\x\. Let u = U\,...,Ut be a list of 
new variables and J be the ideal of fC[u, x\ generated by 1 and h\ — U\, . . . , ht — Ut ■ 
Let G be a Grobner basis of J with respect to the lexicographic order u\ ~< • • • -< 
Ut -< X\ -< • • • -< x n , and {gi, . . . , g v } = Gfl IC[u]. Then: 

i) {gi, . . . , g v } is a Grobner basis of J fl IC[u] with respect to the lexicographic 
order u± -<■••-< u t , 

ii) the fC- isomorphism ui G IC[u\ i — > hi £ JC[h] maps J D IC[u] onto X fl IC[h]. 

In our case, we have I = pK[x, y, z] for some polynomial p, and IC[h] = IC[f, g, s] 
where /, g is a generating system of K,[x,y,z] x and s is a local slice of X. Let 
ui,v,2,U3 be new variables and J be the ideal of K.[ui, u%, 113, x, y, z] generated 
by Pi f — ui, g — 112, s — U3. Let G be a Grobner basis of J with respect to the 
lexicographic order u\ -< 1*2 -< U3 -< x < y ~< z and G\ = G fl /C[iti, U2, U3]. 
According to proposition 3.4, the ideal plC[x, y, z] fl K,[f, g, s] contains a polynomial 
of the form s + a{f, g) if and only if G\ contains a monic polynomial £(u\,U2, U3) 
of degree 1 with respect to 113. In this case the polynomial we are looking for is 
<?(/,<?, s). 

The following algorithm gives the main steps to be performed for computing a 
minimal local slice of a given locally nilpotent derivation in dimension 3. 



Input: A locally nilpotent /C-derivation X of IC[x,y, z] and a generating system 

f,g of K[x,y,z] x . 
Output: A minimal local slice s of X. 

1: Compute a local slice sq of X. 

2: Write X(so) — p" n ■ ■ -p™ r , where the p^s are primes. 
3: s := s . 

4: for i from 1 to r do 

5: for j from 1 to m, do 

6: Let G be a Grobner basis of I(j>i, f — Ui, g — U2, s — 113) with respect to 
the lex-order u\ -< U2 -< U3 -< .t -< y -< z, and let Gi = G n /C[«i, U2, U3]. 

7: if Gi contains a monic polynomial of degree 1 with respect to 113, say 
1*3 + a(u\, U2) then 

8: Write s + a(f, g) = piS\. 

9: s := 8\. 

10: else 
11: Break. 
12: end if 
13: end for 
14: end for 

Algorithm 1: : Minimal local slice algorithm. 
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4. Computing the rank in dimension three 

An irreducible locally nilpotcnt /C-derivation X of JC[x\,. . . ,x n ] is of rank 1 if 
and only if K.[x±, . . . ,x n ] x = K\ n ~ 1 ^ and X has a slice, see [16]. In dimension 3, 
and taking into account theorem 2.3, an irreducible locally nilpotcnt derivation is 
of rank one if and only if Algorithm 1 produces a slice. Therefore, we only need to 
characterize derivations of rank two. 

Theorem 4.1. Let X be an irreducible locally nilpotent derivation ofK.[x,y, z] and 
assume that rank(X) ^ 1. Let us write K,[x,y,z] x = K\f,g\ and S x = cJC[f,g]. 
Then the following are equivalent: 

i) rank(X) = 2, 

ii) c = £{u), where £ is a univariate polynomial and u is a coordinate of K.[f,g], 
Hi) c — t(u), where u is a coordinate of JC[x, y, z] . 

Proof, i) =>• ii) Assume that rank(X) = 2 and let u,v,w be a coordinate system 
such that X(u) = 0. The /C-derivation X is therefore a K, [u] -derivation of /C[u][u, w], 
and since K.[u] is a UFD there exists p £ IC[x, y, z] such that fC[f, g] = K[u,p\. This 
proves that u is a coordinate of lC[f,g]. Let us now view X as /C(u)-derivation 
of K(u)[v,w]. Since X is irreducible and according to theorem 2.2 there exists 
s = h ^( u ™^ such that X(s) = 1, and so X(h) = k(u). Let c be a generator of S x . 
Then c | k(u), and since K,[u] is factorially closed in K\u, v, w] we have c = £{u) for 
some univariate polynomial I. 

ii) =>• Hi) Assume that c = £(u), where u is a coordinate of /C[/, g] and write 
lC[f, g] = K\u,p\. Let s be such that X(s) = c. If we view X as /C(u)-derivation of 
K,(u)[x, y, z] then K.(u)[x, y, z] x = JC{u)[p\ and X(c~ 1 s) = 1. By applying lemma 
2.1 we get K,(u)[x, y, z] = IC(u)[p, s]. From theorem 2.4 we deduce that u is a 
coordinate of /C[x, y, z\. 

Hi) =*> i) Wc have X(c) = £'(u)X(u) = 0, and so X(u) = 0. On the other 
hand, since u is assumed to be a coordinate of IC[x, y, z] we have rank(X) < 2. By 
assumption we have rank(X) ^ 1 and so rank(X) = 2. □ 

The conditions of ii) in theorem 4.1 are in fact algorithmic. Indeed, it is algo- 
rithmically possible to check whether a given polynomial in two variables is a coor- 
dinate, see e.g. [1, 4, 11, 26]. We will use the algorithm given in [11], but it is worth 
mentioning that from the complexity point of view the algorithm given in [26] is 
the most efficient as reported in [22]. On the other hand, condition c = £{u) may be 
checked by using a special case, called uni-multivariate decomposition, of functional 
decomposition of polynomials, see e.g., [18]. It is important to notice here that 
uni-multivariate decomposition is essentially unique. Namely, if c = £(u) = £\(ui), 
where u and u±, are undecomposable, then there exist /i S K.* and v S /C such that 
U\ = f-iu + v. Due to the particular nature of our decomposition problem it seems 
more convenient to use the following proposition. 

Proposition 4.2. Let c(x) £ /C[xi, . . . , x n ] be nonconstant and u = ui,...,u n be 
a list of new variables. Then the following are equivalent: 

i) c fe) = ^(j/ife)); where £ is a univariate polynomial and y\ is a coordinate of 
/C \x\ , . . . , x n ] , 

ii) yi{x) — yi(u) | c(x) — c(u) and y\ is a coordinate of K[x\, . . . ,x n \. 

Proof, i) =>• ii) Let £ be a univariate polynomial and t,t' be new variables. Then 
we have t-t' \ £(t) - £{t'). This shows that y x {x) - yi{u) \ £(y(x)) - £(y(u)). 
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ii) => i) Let i)2, . . . , y n be polynomials such that y = yi, . . . , y„ is a coordinate 
system of IC[x], and let Vi = yi(u). Then v = v\, . . . , v n is a coordinate system of 
K,[u\. Let us write c(x) — c(u) = (yi{x) — yi(u))A(u, x) and c(x) = £(y). Then we 
have 

(4.1) e(y)-£(v) = (yi-v 1 )B(v,y). 

Let us now write t{y) = J2 a a a{yi)%l2 2 ' ' ' Vn" ; where a = (a2, • • . , a n ). After 
substituting i/i to V\ in the relation (4.1) we get 

a a 

Taking into account the fact that t>2, ■ ■ ■ , u n , 2/2, ■ • • , y-n are algebraically independent 
over JC[yi] we get a Q = for any a / 0. This proves that i(y) is a polynomial in 
terms of y\. □ 



Input: A locally nilpotent /C-derivation X of K[x, y, z] and a generating system 

f,g of K.[x,y,z] x . 
Output: The rank of X . 

Write X = a\d x + a2d y + a^d z . Let c\ — gcd(ai, 02, 03) and X = c{y. 
By Algorithm 1, compute a minimal local slice s of y and let c = y(s). 
if c is a unit then 

rank(X) = 1. 
else 

Compute a factorization of c(f,g) — c(ti,t2) in g, ti, t^. 
if no factor of c(f, g) — c(t 1 ,t 2 ) is of the form u(f, g) — u{b\, t 2 ) then 

rank(X) = 3. 
else 

if u is a coordinate of /C[/, g] for a factor of the form u(f, g) — u(ti, ta) °f 
c(/,ff) ~ c(ti,t 2 ) then 

ranfc(Af) = 2. 
else 

rank(X) = 3. 
end if 
end if 
end if 

Algorithm 2: : Rank algorithm. 



As defined in [7], a derivation X of rank r of K.[x] is rigid if for any coordi- 
nate systems y±,...,y n and zi,...,z n such that /C[?/i, . . . , y n -r] £ /C[x]'*' and 
/C[«i, . . . , Zn- r ] C /C[x] we have fC[yi, . . . , y n -r] = fc[zi, ■ ■ ■ , z n -r]- The main 
results of [7] lie behind the fact that locally nilpotent derivations in dimension 3 
are rigid. In fact only the rank two case is nontrivial since in general derivations of 
rank 0, 1 and n are obviously rigid. The characterization ii) of rank two derivations 
given in theorem 4.1 gives a more precise information. Indeed, it tells that if a 
coordinate of K.[x,y,z] belongs to K,[x,y,z] x then it may be found by decompos- 
ing the generator of the plinth ideal S x . The fact that rank two derivations are 
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rigid is then an obvious consequence of the uniqueness property of uni-multivariate 
decomposition. 

In case a given derivation X is of rank 2, Algorithm 2 does not produce a coordinate 
system u, v, w such that X(u) = 0. Computing such a coordinate system, which 
is equivalent to obtaining an algorithmic version of theorem 2.4, is a necessary 
step towards our ultimate goal, namely algorithmically classifying locally nilpotent 
derivations in dimension 3. 

5. Comments on implementation 

Before implementing algorithms for locally nilpotent derivations of K\x,y, z] we 
must first specify how such objects are to be concretely represented. Any chosen 
representation should address the two following problems. 

(1) Recognition problem: Given a /C-derivation X of K.[x,y,z], check whether 
X is locally nilpotent. 

(2) Kernel problem: Given a locally nilpotent /C-derivation X of K,[x, y, z], com- 
pute /, g such that K.[x, y, z] = K.[f, g]. 

Usually, a derivation X of K,[x, y, z] is written as a K. [x, y, z] -linear combination 
of the partial derivatives d x ,d y ,d z . However, with such a representation, the recog- 
nition and kernel problems are nowhere near completely solved. To our knowledge, 
only the weighted homogeneous case of the recognition problem is solved, see [15]. 
One way to go round this hurdle is to opt for another representation. 

The Jacobian representation gives another alternative to represent locally nilpo- 
tent derivations. Indeed, any irreducible locally nilpotent /C-derivation of /C[x, y, z] 
is, up to a nonzero constant in /C, equal to Jac(/, <?, .), see [8]. In order to check 
whether a Jacobian derivation X = Jac(/, g, .) is locally nilpotent it suffices to 
check that X d+1 {x) = X d+1 {y) = X d+1 (z) = 0, where d = deg(/)deg(g), sec [14]. 
However, it is still not clear how such a representation could help in solving the 
kernel problem. Nevertheless, we may always check whether this ring of constants 
is generated over /C by /, g by using van den Essen's kernel algorithm [14] . 

Due to the above discussed issues, we have restricted our implementation to the 
case of derivations of IC[x, y, z] represented in a Jacobian form, say Jac(/, g, .), and 
whose ring of constants is generated by /, g. The computer Algebra system we used 
for implementation is Maple 9. 
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